- 1 - 



Compressible streaming instabilities in rotating thermal viscous 

objects 

A. K. Nekrasov 

Institute of Physics of the Earth, Russian Academy of Sciences, 123995 Moscow, Russia 
e-mail: anatoh . nekr assov@t-onhne . de 

Abstract. We study electromagnetic streaming instabilities in thermal viscous regions of 
rotating astrophysical objects, such as, protostellar and protoplanetary magnetized accretion 
disks, molecular clouds, their cores, and elephant trunks. The obtained results can also be 
applied to any regions of interstellar medium, where different equilibrium velocities between 
charged species can arise. We consider a weakly and highly ionized three-component plasma 
consisting of neutrals and magnetized electrons and ions. The vertical perturbations along 
the background magnetic field are investigated. The effect of perturbation of collisional 
frequencies due to density perturbations of species is taken into account. The growth rates 
of perturbations are found in a wide region of wave number spectrum for media, where the 
thermal pressure is larger than the magnetic pressure. It is shown that in cases of strong 
collisional coupling of neutrals with ions the contribution of the viscosity is negligible. 
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1. Introduction 

In a series of papers, Nekrasov (2007, 2008 a,b, and 2009 a,b), a general theory for 
electromagnetic compressible streaming instabilities in multicomponent rotating magnetized 
objects, such as, accretion disks and molecular clouds has been developed. In equilibrium 
of accretion disks, different background velocities of different species (electrons, ions, dust 
grains, and neutrals) have been found from the momentum equations with taking into ac- 
count anisotropic thermal pressure and collisions of charged species with neutrals. Due to 
velocity differences, compressible streaming instabilities have been shown to arise having 
growth rates much larger than the rotation frequencies. New fast instabilities found in these 
papers have been suggested to be a source of turbulence in accretion disks and molecular 
clouds. 

In papers cited above, the viscosity has not been considered. However, numerical sim- 
ulations of the magnetorotational instability show that this effect can influence on the mag- 
nitude of the saturated amplitudes of perturbations and, correspondingly, on the turbulent 
transport of the angular momentum (e.g., Pessah & Chan 2008; Masada & Sano 2008). The 
viscosity has numerically been shown (Yatou & Toh 2009) can play a crucial role in the per- 
sistence of the long-lived localized clouds observed in interstellar media (Braun & Kanekar 
2005; Stanimirovic & Heiles 2005). Interstellar media, molecular clouds, protostellar and 
protoplanetary accretion disks are weakly ionized objects, where coUisional effects play a 
dominate role. The ratio of the viscosity to the resistivity (the magnetic Prandtl number) 
for astrophysical objects takes a wide range of values. In particular, in accretion disks around 
compact X-ray sources and active galactic nuclei, the magnetic Prandtl number varies by 
several orders of magnitude across the entire disk (Balbus & Henry 2008). Therefore, the 
viscosity is needed to be considered at studying electromagnetic streaming instabilities in 
multicomponent weakly ionized media. 

In the present paper, we explore electromagnetic streaming instabilities in rotating astro- 
physical objects, such as protostellar and protoplanetary magnetized accretion disks, molec- 
ular clouds, their cores, elephant trunks, and so on, taking into account effects of collisions, 
thermal pressure and viscosity. We consider a weakly and highly ionized three-component 
plasma consisting of electrons, ions, and neutrals. The charged species are supposed to be 
magnetized, i.e., their cyclotron frequencies are considered to be much larger than their or- 
biting frequencies and coUisional frequencies with neutrals. We will investigate the vertical 
perturbations along the background magnetic field. The presence of static (in perturbations) 
dust grains can be invoked through the quasineutrality condition. We take into account the 
effect of perturbation of coUisional frequencies due to density perturbations of species, which 
takes place at different background velocities of species. We find expressions for the per- 
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turbed velocity of any species that also contain the perturbed velocity of other species due 
to collisions. For magnetized charged species, we derive the dispersion relation, which is 
solved in the thermal regime when the pressure force dominates the inertia. The conditions 
of strong or weak collisional coupling of neutrals with charged species and the role of the vis- 
cosity will be considered. The growth rates due to different azimuthal velocities of electrons 
and ions will be found. 

The paper is organized as follows. In Section 2 the basic equations are given. In 
Section 3 we shortly discuss the equilibrium state. Solutions for the perturbed velocities of 
species for the vertical perturbations arc obtained in Section 4. The dispersion relation in 
the general form is derived in Section 5. In Section 6 this dispersion relation is solved in 
the thermal regime in the specific cases and unstable solutions are found. In Section 7 we 
give an expression needed for determining of polarization of perturbations. Some problems 
concerning the contribution of the resistivity in the standard MHD are discussed in Section 8. 
Discussion of the obtained results and their applicability to protostellar and protoplanctary 
disks are given in Section 9. The main points of the paper are summarized in Section 10. 



2. Basic equations 

We will consider weakly ionized rotating objects consisting of electrons, ions, and neu- 
trals. Here, we do not treat the presence of dust grains. However, the latter may be involved 
as static (in perturbations) species through the condition of quasineutrality. The electrons 
and ions are supposed to be magnetized, i.e., their cyclotron frequencies are larger than 
their collisional frequencies with neutrals. Self-gravity is not included. We will study one- 
dimensional perturbations along the background magnetic field Bq. Then the momentum 
equations for species in the inertial (nonrotating) reference frame accounting for the viscosity 
(Braginskii 1965) take the form, 

^+v,--Vv,- = _vC/-^ + ^E+^v,xB-i.,.„(v,-v„) (1) 

d^Wj X b 6 /ij c?^vj_L 4 



+ K , ,2 I2 + Q^i" 



+ V„ • VV„ = -VU - —— - Vr^j (v„ - V,) + -— — , (2) 



at mnUn ^-^ \ oz^ 3 az^ 

3 

where the index j — e,i denotes the electrons and ions, respectively, and the index n denotes 
the neutrals. In Equations (1) and (2) qj and the charge and mass of species j 

and neutrals, Vj,„ is the hydrodynamic velocity, rij^n is the number density, Pj^n — nj^nTj,n 
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is the thermal pressure, is the temperature, i/jn — ^jrjnnnn {^nj — Ijnfnjnj) is the 
collisional frequency of species j (neutrals) with neutrals (species j). The indices _L and z 
denote directions across and along the magnetic field, respectively. The value is 7j>i =< 
av >jn I {rrij + mn), where < av >j„is the rate coefficient for momentum transfer, and fij n = 
"^Tj nl ^jn,nn are the coefficients of the kinematic viscosity (z/^n is the neutral-neutral collisional 
frequency), where VTj,n = ^j,n/^j,nY^'^ is the thermal velocity Further, ujcj = ^jBo/rrijC is 
the cyclotron frequency and Tjn = vj^- Other notations are the following: U — —GM/R 
is the gravitational potential of the central object having mass M (when it presents) , G is 
the gravitational constant, R — {r^ + z'^Y^'^, E and B are the electric and magnetic fields, 
h — 'B/B, and c is the speed of fight in vacuum. The magnetic field B includes the external 
magnetic field Boerrt of the central object and/or interstellar medium, the magnetic field 
Bocur of the stationary current in a steady state, and the perturbed magnetic field. We 
use the cylindrical coordinate system {r,9,z), where r is the distance from the symmetry 
axis z, and 6 is the azimuthal direction. We assume that the background magnetic field is 
directed along the z axis, Bq = ^ozext + ^ozcur- In Equation (1), the condition Ucj ^ I'jn is 
satisfied for the viscous terms. For unmagnctized charged particles of species j, Ucj ^ i^jn, 
the viscosity coefficient has the same form as that for neutrals. We adopt the adiabatic 
model for the temperature evolution when Pj^n ~ 't-J^, where 7a is the adiabatic constant. 

The other basic equations are the continuity equation, 

-^ + V-n,-„v,-„ = 0, (3) 

Faraday's law, 
and Ampere's law, 

V X B (5) 
c 

where j — We consider the wave processes with typical time scales much larger 

than the time the fight spends to cover the wavelength of perturbations. In this case, one 
can neglect the displacement current in Equation (5) that results in quasineutrality both for 
the electromagnetic and purely electrostatic perturbations. 



3. Equilibrium 



We suppose that electrons, ions, and neutrals rotate in the azimuthal direction of the 
astrophysical object (accretion disk, molecular cloud, its cores, elephant trunk, and so on) 
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with different, in general, velocities Vj^no- The stationary dynamics of light charged species, 
electrons and ions, is undergone by the effect of the background magnetic field and collisions 
with neutrals. In their turn, the neutrals also experience the collisional coupling with charged 
species influencing on their equilibrium velocity. Some specific cases of equilibrium have 
been investigated in papers by Nekrasov (2007, 2008 a,b, 2009 b), where the expressions for 
stationary velocities of species in the gravitational field of the central mass have been found 
at the absence of coUisions as well as with taking into account collisions for cases of weak 
and strong collisional coupling of neutrals with ions. 

Due to different stationary velocities of charged species, the electric currents are gener- 
ated in the equilibrium state. 



4. Linear regime 

In the present paper, we do not treat perturbations connected with the background 
pressure gradients. Thus, we exclude the drift and internal gravity waves from our consid- 
eration. Wc take into account the induced reaction of neutrals on the perturbed motion 
of charged species. The neutrals can be involved in electromagnetic perturbations, if the 
ionization degree of medium is sufficiently high. We also include the effect of perturbation of 
the collisional frequencies due to density perturbations of charged species and neutrals. This 
effect emerges when there are different background velocities of species. Then the momentum 
equations (1) and (2) in the linear approximation take the form. 



"""1 ^ , N 9''vj, X bp 6 fij 9V,ai 4 SV 

■—-,n (V.0 - v„o) + ^^^^ + 5 ^^i^ + 3^^-^ 



jzl 



^nl 2 ^^nl / \ '^jl 0* / \ 

~W ^ '""'^^ 2^i^n,(Vnl-V,i)-_^— I/„^.(v„o-V,o) (7) 

J 3 ■' 

f (9Vni 1 d\nzl 

[ dz^ ^3 dz^ 

where Csj^n — {laTj,no/fnj,n)^^'^ is the sound velocity (7a = 3 in the one-dimensional case), 
^'jn = Iju^nnnOj ^nj = 7jn^i^jo- The tcrms proportional to Uj* and z/,°* describe the 
effect of perturbation of the collisional frequencies due to density perturbations. The index 
1 denotes the quantities of the first order of magnitude. The neutrals participate in the 
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electromagnetic dynamics due to coUisional coupling with the charged species, mostly with 
ions (in the absence of dust grains). However, we keep, for generality, collisions of neutrals 
with electrons in Equation (7). 

The continuity Equation (3) in the linear regime is the following: 

%^ + n,,„o%^ = 0. (8) 

We further apply the Fourier transform to Equations (6)- (8), supposing the perturba- 
tions of the form exp(ikzZ — iut). Using Equation (8) and solving Equation (7), wc find 
the expressions for components of the Fourier amplitude Vnik, where (and below) the index 
k = {A;^,a;}, 

— iLOn±Vnlrk — ^njVjlrk, (9) 



-itiJn±Vniek 



EV^ * kz{VnO-Vjo) 
l^njVjWk - Z^^nj '"nzk, 



UJ 

J 3 



—ioJnzVfilzk — J^njVjlzk- 



3 

Here and below, the index by v'^^ and is, for simplicity, omitted. In Equations (9), the 
notations are introduced, 

3 

a;„_L = uJ + ll'n + ^^Xnz, 

<-^nz = VlVn + ^Xnz, 

UJ 

where i/„ = Y^- u^j, Xnjz = llJ-njkl- 

Now we substitute the expressions for v^i^ and rij^nik = n-j^nokzVj^nizk/^ in Equation 
(6). Then we obtain the following expressions for components of v^i^: 



-iuJj±VjXrk = —Eirk+(^jVjWk + i-^y^,T^nlVllrk, (10) 

rrij i^n±^ 

—ZUJj±Vji0k — Ei0k — UJjVjirk + I / ^J^nlVlWk — i / .PjlzVllzk, 

'rrij ^ ^ 

ilzk — I ^Izk + 'r^z ^lek ] + ^ / , ^nlVllzk- 



-lUjzVji 
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Here the notations are introduced: 



ujjj_ = a; + iujn + iXj±2, 

UJj^ ^ U ^ + iVjn + iXjz, 

UJ 

ujj — ujcj — 

Pjlz - ^i'^^n/-- — + — 

where Xjj_i = Hjk'^/ucjTjn, Xj±2 = ^■'^l^jk^/'-^cj'^jn- ^^^^ Equation (10), we have 

substituted Bi^k for -E'lefc; using Faraday's equation (4), Bi^k — —UzEiq^ {Biek — f^zEirki 
Bizk — 0), where Uz — kzc/uj. 

From Equation (10) for Vjizk, we can find the longitudinal (along Bq) velocities of 
electrons and ions, 

DzVjizk = irizhjzEiek + iajzEuk- (11) 



Here, 



Dz — ^ez^iz^nz + ^iz^en^ne + ^ 

- A _ A 

O'ez — ^iz ^en^ni ; 

me mi 

rUe c rui c 

rrii rrie 

Oiz — Oiez ^in^ne i 

rrii c me c 



Let US find now the transverse (across -Bq) velocities of charged species. When solving 
Equations (10), we will consider the case in which the charged species are magnetized, i.e., 
we will suppose that the following conditions are satisfied: 

Ucj > UJj±, iyj„Un/uJn±, Xj±l- (12) 

Conditions (12) signify that the Lorentz force is dominant. Then, using Equation (11), we 
find solutions, 

Qj 1 / Qj Tlz \ ^jz 

Vjirk — —i n-bjEirk H — aj + — Eigk H —Ei^k, (13) 

rrijUJ^j ujcj \mj Dz ' J cUcjDz 

VjWk — CljEirk — i o-bjEigk- 
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Here, 





UJcj UJn± 




Qe Qi 

^jez ~r ^iiz-i 
Trie TUi 




Qe VeO qi Vio 
Ijez + 7ji2 ) 

me c rrii c 


Ifjez — 


Pjez^^iz Pjiz^in^nei 


Ijiz — 


Pjiz^ez Pjez^en^ni- 



In the next section, we calculate the perturbed electric current and obtain the dispersion 
relation. 



5. Dispersion relation 

Prom Equations (4) and (5) we obtain. 



Jlzk 



UJ 



-Ji±ifc, 



(14) 



0. 



Using solutions (11) and (13), we can calculate the electric current ji^ = J2j 1j''^jo^jik + 
X^^- g^njifeVjo- Substituting this current in Equations (14), we will find the following set of 
equations for determining the components of the perturbed electric field Ei^ = {Eigk, Eirk, Eizk) 



AfeEifc — 0, 



(15) 



where the matrix is equal to 



-i (£2 + £x-L + ^A^) , 
£4, 



I {€2 + £x±) 

nl - £1, 



£& 

-i£5z 
-£3 



The components of the matrix are the following: 



^1 - Z.^^^'^2--^,^3--^2.'^p.— «^-^4-^Z.'^p.— ^^^^^ 
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£5 — — pr / .^pj Ojz,£Q — —fz- / .^pi «i2) 

UjDz ^ Qj C CuDz ^ Qj C 

_ 1 u^^.m ^ _ nz ^^pjm 

^ UJUJ^j UUz ^ UJcj Qj UJVz UJcj Qj 

where cOpj = (^inrijoq'^ / nij) is the plasma frequency and the index d denotes dust grains 
(if they present and are static in perturbations) . 

The dispersion relation is obtained by equating the determinant of the matrix Aj. to 
zero. As a result, we obtain, 

[{nl - £1) £3 + £365 + e^se] {nl - £1) - (£2 + e^±) [(£2 + e^±) £3 + £3£az + SiEsz] = 0. (16) 

In this equation, it is easy to see that the value £2+£^x-L general, small in comparison 

to the value — £1. Therefore, we consider below the case in which the following conditions 
are satisfied: 

{nl-e,) » (£2 + £x±)> (17) 

(£3£5 + £4£6) {nl - Si) > (£2 + S^_l) {SsExz + S^SSz) ■ 

Then Equation (16) will take the form, 

{nl - £i) £3 + £3£5 + £4£6 = 0. (18) 

We do not consider the damping Alfven perturbations — £i = 0. When Equation (18) is 
satisfied, the component of the electric field Ei^^zk ^^'^ ^irk Eiq^. The second and 
third equations of a set (15) determine the components Eirk and Eizk, 

{nl - £i) £3-Ei^fe = i [is2 + ^3 + essxz + e^esz] Eieu, (19) 
^sEizk — e^Eiek- 

Below, we will investigate solutions of Equation (18). 

6. 6. Dispersion relation (18) and its solutions in the specific case 

We further will not take into account the presence of dust grains. Then the quasineu- 
trality condition has the form ge^eo + ?i^io — 0. Using this condition, we will calculate the 
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quantities £i, £3, £4, £5, and Sq. We will be interested in the case in which a; <S + Xnz- 
When Xnz < ^'n there is strong coUisional coupling of neutrals with charged species in the 
transverse direction to the z axis. This coupling along the magnetic field depends on the 
relation between A;^c^„ and uju^ [see Equations (9)]. If Xnz ^ '^n, then the neutrals have 
a weak coupling with charged species. We suppose also that uj ^ Ujn, Vnn (this inequality 
is wittingly satisfied in weakly ionized plasma if Xnz ^ and u) -C Un). In this case, the 
contribution of the viscosity X3,nz in the expressions for ujj^nz is much smaller in comparison 
to the thermal term klc^j j^/u; and can be ignored. We further will investigate the case in 
which k1(?gj ,^ ^ uj"^ when the thermal pressure dominates the inertia. Below, we will obtain 
the dispersion relations for two regions of small and large wave number kz- 



6.1. Dispersion relation (18) for small k^ 

At first, suppose that the following condition is satisfied: 



where 



ujd^ » klc%cl^, (20) 



Under conditions at hand we will find. 



3 



2 2 ' 



e I ; 



- ^le^zk l 



■ pe z'^z J 

£4 = -^ — jj^di, 



£5 — — * — f; "2 H F;~^t'2 

uDz OJ UjJJz c 



£6 = £4 + -^^U,yi , (21) 



where 



,2 



ui — Cbi (Ze, 02 — —^Cli T"e) 

C QelTli C Qe'nT'i & 
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_ qinie 

_ ^ViO qirrie Vio 

c q^nii c 

In the expression for £i, we have supposed that ^ Xnz that is compatible with condition 
(20). 

The value Dz has the form, 

Dz = g^, (22) 



where 



2 _|_ 2 _|_ 2 



Substitute now the expressions (21) and (22) in Equation (18). Then we obtain the 
following dispersion relation under conditions given above: 

inl - Si) { —di - -de ) g + [ikldedi + U {Ven^nidi + Vin^^nede)] (^^iO " V^of = 0. 

(23) 

6.1.1. Solution of dispersion relation (23) for sufficiently leirge kz Let us find a 
solution of Equation (23) in the case 

kldedi > ixl {l^en^nidi + l^in^nede) ■ (24) 

Then Equation (23) takes the form 
where 

2 2 dgdi 



J yrue ' mi '^J 



The sign | | denotes an absolute value. We see that due to the last term on the left-hand 
side of Equation (25) there is an unstable solution. 
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The solution of Equation (25) in the region o;^ <^ k^c^ is equal to 

7 = /i2 



2K0 - Viof 



(26) 



where 7 =Im iv is the growth rate. Thus, for kz such that kz > k* = h"^ (feo — Vio)'^ /cac^, the 
growth rate is maximal and independent from the wave number. In the region a;^ ^ ^1'^% 
we obtain, 



7 



{VeO - Vio) 



2/3 



(27) 



This growth rate increases with increasing of in the region k^ < k*. 



6.1.2. Solution of dispersion relation (23) for sufficiently small k^ In the case 



(2^ 



Equation (23) has a solution. 



UJ — 



,2 2 ('^e0 - VioY 



1/2 



(29) 



where 



2 2 

S = UJ 



pe 



Y me ' rrii ^ J ^ 



We see that perturbations with the wave number kz < s\veQ — Vio\/c will be unstable. 



6.2. Solutions of dispersion relation (18) for large k^ 

Now we consider the dispersion relation (18) in the case of large k^ when 

.22 2 



(30) 



Then we obtain the following expressions: 

,2 



UJ, 



£3 



^3J^ -z-sn \ -St 



qiiri 



e „2 



-13- 



Si 



qeUli c 



£5 



^le^l 4 2 ^ _ 
3£) '^z'-sn \ ^2 s« 



e0^2 2 \ , ^^^^e^ ^iO - ^eO 



£6 — £4 H ) 



7-» '"2 2 2 2 



a; 



3 se St sn' 



(31) 



The expression for £1 in the case Un ^ has the form (21). In the opposite case, z/„ -C Xnz, 
we obtain, 



(32) 



Substituting expressions (31) and Ei from (21) in Equation (18), we obtain in the case 
^ Xnz a dispersion relation which has a solution. 



UJ — 



,2 _ (^eO -Viof 



1/2 



(33) 



where Cg = [3 (lo'el^'i + QiTe) /\qe\mi]^^'^ is the ion sound velocity. Using expressions (31) and 
(32), we find in the case Xnz ^ J^n the following solution of the dispersion relation (18): 



UJ = I 



Yl ~f ^^3n + X3I-2) 



-1 



UJ, 



— pi I \2 T 2 2 

(^;eo - Vio) - k^c 



(34) 



The perturbations (33) and (34) with wave numbers such that k^ < kth, where kfh — 



are unstable due to different equilibrium velocities of electrons and ions. 



7. On calculation of Eij.^ 

To find the radial component of the electric field Ei^k defined by Equation (19), it is 
necessary to calculate the value e^ex^ + e^esz- The terms e^z and Sx^. have been arisen due 
to perturbation of coUisional frequencies proportional to the number density perturbations. 
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Using the expressions for these terms containing in Equation (15) and the expressions for £3 
and £4 from Equations (20), we obtain, 

£2.£\z + £4^5^ = I I > ^ 

uj^Dl mi\ c ) ^ Ucj Qj 

where Pjiz, I = e,i, is defined in Section 4. We see that the value EsExz + SiEsz ~ ("t^eo — 
Vio){veo — f„o) + {veo — fio)('i^jo ~ "^no)- Thus, polarizatiou of perturbations also depends on 
the difference of equilibrium velocities of species. 



8. On the resistivity in the standeird MHD 

In the papers by Nckrasov (2007, 2008 a,b, and 2009 a,b) and in the present paper, we 
study streaming instabilities of multicomponent rotating magnetized objects, using the equa- 
tions of motion and continuity for each species. From Faraday's and Ampere's laws wc obtain 
equations for the electric field components (see Equation (15)). Such an approach allows us 
to follow the movement of each species separately and obtain rigorous conditions of consid- 
eration and physical consequences in specific cases (see, e.g.. Section (9)). This approach 
permits us to include various species of ions and dust grains having different charges and 
masses. In some cases, the standard methods used in the magnetohydrodynamics (MHD) 
leads to conclusions that arc different from those obtained by the method using the elec- 
tric field of perturbations. Indeed, let us consider the magnetic induction equation. For 
simplicity, wc take a two-component electron-ion plasma embedded in the background mag- 
netic field directed along the z axis. The magnetic induction equation with the resistivity is 
obtained from the momentum equation for the electrons at neglecting the electron inertia, 

= -— E —VeXB-Uei{Ve-Vi), (35) 

rUe rrieC 

where Uei {i^ie) is the electron-ion (ion-electron) coUisional frequency, and — e is the electron 
charge. Replacing Vj — Vg by the current j/en (rig = rii = n), applying Vx to Equation (35), 
and using Equations (4) and (5), we obtain the well-known magnetic induction equation, 

= V X Ve X B+T^^V^B, (36) 

where rjm = ^eiC^ /^pe is the coefficient of the magnetic diffusion or the resistivity. The 
momentum equation for ions has the form, 

e „ e ^ , , , 

— ^ = — E+ Vi XB-Uie Vi - Ve , 37 

at rrii rriiC 
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where d/dt — d/dt+Vi-Vvi. It is seen from Equations (35) and (37) that if the electrons and 
ions are magnetized, then Ve_L ~ Vj^ (the sign _L denotes the transverse direction relatively 
to the magnetic field). Usually, the electron velocity in Equation (36) is substituted by the 
ion (or neutral) velocity. As a result, one solves the standard two MHD equations, 

= -LvxBxB, (38) 
dt An ^ ' 

— = V X V X B+77^V'B, 

where p — rriin, v = Vj. 

Prom Equations (38) in the linear approximation, one obtains the dispersion relations, 

+ iujTjmk'^ - k^c\ = 0, (39) 

for magnetosonic waves, and 

oj^ + iojrjmk^ - klc\ = 0, (40) 

1 /9 

for Alfven waves. Here, k'^ — k^^- k^^ ca — {Bl/Anp) is the Alfven velocity. 
Two consequences are followed from Equations (39) and (40): 

1) . Magnetosonic and Alfven waves are damped due to the resistivity. 

2) . The resistivity is isotropic since it is proportional to the total wave number. 

We derive now the dispersion relation, calculating the perturbed velocities of electrons 
and ions through the electric field. For magnetized electrons and ions, we find from Equations 
(35) and (37), 

enc ijj ( ijj \ 

jlrk = -B ( -iEirk H E-^gk , (41) 



enc 


UJ 


Bo 


^ci 


enc 


UJ 


Bo 






e 



jwk — -f^ ( —iEidk Eirk 

Velz = Eu,Viiz^O. 

niel^ei 

We see that the coUisional frequency is absent in the transverse current and determines 
only the electron velocity along the background magnetic field. Using Equations (4) and 
(5) and neglecting the small terms proportional to uj/iOci in the round brackets of Equations 
(41), we obtain the following dispersion relation: 

(a;^ - k'^c\) (o;^ + iur^mkl - klc\) = 0. (42) 
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We can conclude from Equation (42): 

1) . Magnetosonic waves are not damped because they have no field 

2) . Alfven waves arc damped due to resistivity. The resistivity is proportional to k'j_ 
and is anisotropic. When k± = Alfven waves also are not damped. 

These results differ, in principle, from the one that follows from Equations (39) and 
(40). From our viewpoint. Equation (42) takes into account the physical mechanism of the 
coUisional damping correctly. 

9. Discussion 

The neutrals participate in the electromagnetic perturbations only due to collisions with 
the charged particles, electrons, ions, and dust grains. One can say that the neutrals are as 
passive agent, ballast making difficult perturbations of charged species. At the same time, 
the latter are only active agents generating electromagnetic perturbations. Therefore, an 
adequate description of multicomponent plasmas including the neutrals is to express the 
neutral dynamics through the dynamics of charged species and substitute induced velocities 
of neutrals into coUisional terms of the momentum equations for charged species. Then using 
Faraday's and Ampere's laws, we can derive the dispersion relation in the linear approxima- 
tion and/or investigate nonlinear structures. 

From expressions for induced velocities of neutrals, one can easily to find rigorous condi- 
tions when the neutral dynamics is important or not, i.e., when there is strong (or sufficiently 
strong) or weak coUisional coupling of neutrals with charged species. In the last case, the 
neutrals are immobile in the electromagnetic perturbations. 

Let us write out conditions, which have been used for obtaining the dispersion re- 
lation (18), and consider parameters of astrophysical objects, for which these conditions 
can be satisfied. We will apply our results to protostellar and protoplanetary disks. We 
have assumed that the electrons as well as ions are magnetized, oOcj 2> Pjn- The condition 
u}ce{> 0) ^ Uen is, in general, satisfied in astrophysical objects (Wardle & Ng 1999). The 
rate coefficient for momentum transfer by elastic scattering of electrons with neutrals is 
< au >e„= 4.5 X 10-9(Te/30 K)^/^ ^^^3 g-i p^aine et al. 1983). As for ions, we should 
consider parameters to satisfy conditions Ud ^ I'in and 3> c^j ^ c\ or c^j ^ (the 
last two conditions see below). Thus, the magnitude of the magnetic field must be in some 
limits. We take the standard values = SOnip and m„ = 2.33mp {nip is the proton mass). 
The rate coefficient for momentum transfer < au >inis equal to < au >in— 1.9 x 10~^ cm^ 
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S'"^ (Draine et al. 1983). Then, for example, from conditions cuci ^ J^in and c^j ^ c\ we 
obtain (3n„Tj)^^^ 3> -Bq ^ 0.43 x 10~^^n„ (Tj is in the energetic units, Bq is in G, and n„ is 
in cm^'^), where we have used qi = —Qe- Note that it is followed from the last inequalities 
that the number density of neutrals is limited for a given ion temperature. If we take Tj 
(K) = 700 K, then wc obtain n„ -C 1.57 x 10^^ cm~^. In this case, the neutral mass density 
Pn = mnTin "C 6.1 X 10~^^ g cm~^. This condition is applicable for early stage of protoplan- 
etary disks or in the surface layers, where the density is lower and temperature is higher. 
In the dense inner parts of a disk, where p„ ~ 10"^'^ — 10~^ g cm~^ (Hayashi et al. 1985; 
Wardle & Ng 1999), the ions can be unmagnetized and their viscosity will be of the same 
form as that for neutrals. Our model is not applicable to such dense regions. 

For rotating protostellar cores of molecular clouds (protostellar disks) , we will adopt the 
following parameters: n„ = 10^ — 10^ cm~^, ni/rin — 10~^ — 10~^ (e.g., CaseUi et al. 1998; 
Ruffle et al. 1998; Pudritz 2002), and Bq = 10 //G (e.g., Goodman et al. 1993; Crutcher et 
al. 1999; Caselli ct al. 2002). For this magnetic field we obtain ujce — 1.76 x 10^ s~^ and 
Uci = 3.19 X 10-3 s-^ 

Now, we give some relationships which are useful at analysis of conditions of consid- 
eration. For protoplanetary disks, we take Tg = 700 K. Then we obtain Ven/i^in — 158.73, 

Vnel^en = 2.34 X 10~'^ne/nn, I'ni/l^in = 12.88ni/n„, Uni/^ne = 3.47 X lO'^ni/Ue, I^en/ = 

12.32n„/ni, rriiUin/mei'en = 3.47 x 10^, Unn/'^in = 6.94 (< au >„„~< au >j„). For pro- 
tostellar disks, wc take Tg = 70 K. Then Ven/^in = 50.19, Vni/^ne = 1-1 x lO^nj/ng, 
T-'en/^n = ^-^nn/rii, and iTLiVin/ iTLf.Vf.n = 1.1 X 10'^. Wc SCO that duc to more low electron 
temperature, the role of electron-neutral collisions in comparison to the ion-neutral colli- 
sions becomes less in protostellar disks than that in protoplanetary disks. It can be shown 
that inequalities (17) are wittingly satisfied under conditions used in Section 6. 

For small kz (Section 6.1), we will suppose that the ionization degree is ni/rin <^ 
2.88(0.91) X lO^'^ for protoplanetary (protosteUar) disks (for the large kz (Section 6.2) this 
condition is not necessary (see below)). In astrophysical objects, such as interstellar medium, 
molecular clouds, and protostellar and protoplanetary disks the ratio satisfies in gen- 

eral this condition. In this case, the value — c^sn^en (see inequahty (20)). The value di has 
the similar form, d; = c?g^Vj_m blI rii/nn <^ 1. We adopt that Tf. ^ Ti ^ = T. Then, for the 
thermal effects to be important for all species, k^clj^^ ^ a;^, the condition /c^cfj 3> cu^ will be 
sufficient. The condition (20) for ions and electrons will be satisfied, if ujUen ^ ^zC^e- The 
condition (24) takes the form A;^c^„ 3> uiUni- In this case, there is weak collisional coupling 
of neutrals with ions along the magnetic field (see expressions (9)). It is followed from the 
last two inequalities that the condition mgi'en ^ nT'n^ni must be satisfied for solutions (26) 
and (27) to be realized. It was supposed above when we have obtained de — c^n^en- Under 
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conditions at hand, the value g — c^snyin^en and — cUpJuin (see Equation (25)). Then 
solution (26) takes the form, 

^{j^^o^::^^ (43) 

Collecting all conditions given above, we find that solution (43) is satisfied for wave numbers 
in the region 

^^^^ ^ t.2 ^ ^^"^ 1 

Note that we consider the case in which c^^ » c\. Using condition ^/vni <S 1, we obtain 
from (43) that inequality 16.29 jt'eo — fiol ^ ('^jo is in cm s~^) must be satisfied. We 

see that for protostellar disks this condition is not realized. For protoplanetary disks we 
take n„ = 5 X 10"*^^ cm~^. Then we obtain |feo ~ "^iol 4.34 x 10'* cm s~^. We further take 
\veo — Viol = 3 X 10^ m s^-*^, rii/nn = 10^^, and Bq = 0.25 G. In this case, Upi = 5.37 x 10'^ 
s-\ = 79.75 s-\ = 68.47 s-\ Vni = 8.82 x 10"^ s-\ v^n = 1-09 x 10^ s^^ The Alfvcn 
velocity ca = 5.05 x 10^ m s"^ Csn = 2.73 km s"^ Csi = 7.6 x 10^ m s"^ and Cse = 1-78 x 10^ 
km s~^. The growth rate 7 is equal to 7 = 4.2 x 10~^ s~^. The wave number is in the 
limits, 3.8 x 10~^ m~^ ^ 0.83 x 10~^ m~^. The wavelength of unstable perturbations 
Xz — 27: /kz has a range 

7.57 X 10^ km > > 1.65 x lO^km 

Solution (27) has the form 



7 



2/3 



(44) 



This solution is satisfied for wave numbers in the region 

^3/4 / X 3- 



mm 



{klcs 



\3/2 



> ^ > max 



kz Csi 



kz Csn 



3/2" 



Prom these inequalities we see, in particular, that the 



must be satisfied. We can also conclude that 3> c] 



where A;* = i^mZ^m/^ci l^'eo - t'iol- 

condition \VeO - Vio\ /Cgi < Vin/iOd muou oa-uionc.^. vvc ..0,11 0,10^. ..^.ii..iuvic uiicn, u^j ^ 

in this case. The neutrals are strongly (weakly) coupled with the ions across (along) the 
magnetic field in perturbations (43) and (44). The viscosity of neutrals along (across) the 
magnetic field is negligible in comparison to the thermal pressure (z/„i). The viscosity of 
electrons and ions is negligible because /c^p^ ^ 1, where pj — VTj/oJcj is the Larmor radius. 

Let us consider solution (44) for specific parameters in protostellar and protoplanetary 
disks. For protostellar disks we take Un — 10^ cm"^ and rii/nn — 10 
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Then we obtain = 
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1.37 X 10~^ and Vm/oJci — 4.29 x 10"^. In this case, the condition \veQ — Vio\ <^ CgiVin/oJci 
takes the form |veo — t'iol 1-03 x 10^ cm s~^, where Cgi = 2.4 x 10^ m s~^. In real situation, 
the value \veo — Vio\ is most probably larger and does not satisfied the last condition. Thus, 
solution (44) does not realized in protostellar disks under conditions at hand. 

Let us determine parameters of protoplanetary disks for which solution (44) is satisfied. 
This solution exists if \veo — Vio \ /cgi < Vin/^d- In the case when collisions of neutrals with 
ions do not influence on the equilibrium velocity of neutrals f„o, the difference |feo — fjol can 
be estimated as \veo — fjol ~ '^no^^in/^ci (Nekrasov 2009 b). Thus, the last inequality takes 
the form, Vno/csi < uOci/^in^ where Cgi = 7.59 x 10^ m at T; = 700 K. We consider the 
region of rotating (pre)protoplanctary disk where Vno ~ 10 km s^^. In this case, we obtain 
condition Uin/(^ci < 7.59 x 10~^. We take Uin/uJci = 3 x 10~^ and Bq = 0.5 x 10"^ G. It is 
followed from this that cUd = 1.6 s~^ = 4.8 x 10"^ s"-^, and n„ = 3.5 x 10^ cm~^. The 
ionization degree we take to be rii/nn — 10~^. Then i/„j = 6.17 x 10~^ s~^. Other parameters 
are equal: LUpi = 4.49 x 10^ s~^, \veo — Vio] = 9 m s~^, = 3.82 x 10^ m s~^. Then we obtain 
kl — 2.05 X 10~^^ m~^. The wave number kz satisfies conditions 



Under parameters at hand, we have 0.48 ^ kz/k* 0.06. Thus, the wavelength of pertur- 
bations is in the band 

5.13 X 10^ km > > 6.38iea;trml0^ km. 

For kz/kl = 0.4 the growth rate (43) is equal to 7 = QM^ni = 3.33 x 10"^ s'^. 

Perturbations with more small kz, k1(?g^ <C uJVni (see inequality (28)), have the growth 
rate (29). When additionally kz <^ s\veo — Vio\/c, where s — {ujpi/csn) {'i^ni/'i^inf^'^ ^ this growth 




rate is 



(45) 



7 = i/, 



m 



C. 



■sn 



The wave number kz is in the region 




It is followed from here that conditions Vni ^ 12. 887 and c^j ^ c\ must be satisfied. In 
these perturbations, the neutrals arc strongly coupled with ions both along and across the 
magnetic field. The viscosity of species is also negligible in the case (45). For parameters 
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given above condition ^ 12. 887 results in a large magnetization, Vm/oJci ^ 10~^, for 
which the value \veo — Vio\ is too small. Thus, solution (45) is not available for protostellar 
and protoplanetary disks. 

Let us now consider the large kz when k^cl^ ^ ujUin (see inequality (30)). Under this 
condition the inequality k^c^^cl^ ^ code is satisfied at rii < n„. In this case, the neutrals 
are weakly coupled with ions in the direction along the magnetic field (see Equations (9)). 
Solutions considered below exist in media where c^- ^ or 127rnjTj ^ Bq. This condition 
is not satisfied for typical parameters in protostellar disks given above. However, near the 
central star, an innermost part of the protoplanetary disk can be highly ionized due to high 
temperature and have a high density. Note that a highly ionized dense plasma disks exist 
around neutron stars, stellar black holes, active galactic nuclei, and white dwarfs (see, e.g., 
Jin 1996 and references therein). The last condition together with cua ^ J^in results in 
(UnriiTi)^^'^ > -Bo > 0.43 x lO-^^n^. We see from here that at T ~ 10^ K and rii ~ Un the 
number density must be Un <^ 2.62 x 10^^ cm~^. If we take n„ ~ 10-*^^ cm~^, we obtain the 
range of the magnetic field, 22.81 Bq 3> 4.3 G. Let, at first, Xnz -C i^ni or /c^c^^ -C I'ni^nn 
when the neutral viscosity does not play a role. Then the neutrals are strongly coupled 
with ions across the magnetic field. Under conditions at hand we see that z/„j ^ 1.86c<j 
should be satisfied. In this case, the solution for uj has the form (33). For kz such that 
kz -C kth = ^pi \veo — t'iol /csC, the growth rate is equal to 



The wave number kz is in the region 



Using parameters given above and taking = 7 G, we find: ujpi = 7.59 x 10^ s ^, 
Uci = 2.23 X 10^ s-\ = 1.37 x 10^ s'^, Vni = 1-76 x 10^ s-\ ca = 3.16 km s^^ 
cai = 0.88 km s^-*^, = 4.06 km s^^, Cgi = 2.87 km s^-*^, and Cgn = 10.3 km s^^. From 
condition z/„j ^ 1.86c<; we can find the limit on the value \veo — Vio\ <^ 4.81 km s~^. In 
the case |feo — fiol = 1 km s~^, we obtain kth = 0.62 m~^.The growth rate (46) is equal 
to 7 = 1.96 x 10^ s~^. The wavelength of unstable perturbations has a range 1.1 x 10^ 
m» Xz ^ 10.06 m. The viscosity of electrons and ions is neghgible since k^pj <^ 1. 

The strong coupling of neutrals with ions in the transverse direction is maintained at 
Xnz < ^ni- When Xnz ^ ^ni or /c^c^^ ^ i^ni^nn the ueutral-iou coupling is weak. In this case. 
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the solution (34) at kz <S kth can be written in the form 
The condition Xnz ^ 7 is satisfied. The wave number is in the region 



From these inequahties we see that c^j ^ c^- (we assume k'^pf < 1). The growth rate (46) 
can be larger than z/„j. At the same time, condition of consideration is 7 <^ z/j„ (see Section 
(6)). Thus, the condition \veo — Vio \ /cg -C I'in/^d must be satisfied. Note that solution (47) 
is the only case when the viscosity of neutrals and ions is important. However, the neutrals 
are weakly coupled with ions in this case and are immobile in perturbations. 



10. Conclusion 

In the present paper, we have studied electromagnetic streaming instabilities in thermal 
viscous regions of rotating astrophysical objects, such as, protostellar and protoplanetary 
magnetized accretion disks, molecular clouds, their cores, and elephant trunks. However, 
the obtained results can be applied to any regions of interstellar medium, where different 
background velocities between electrons and ions can arise. 

We have considered a weakly and highly ionized three-component plasma consisting 
of electrons, ions, and neutrals. The cyclotron frequencies of charged species have been 
supposed to be much larger than their coUisional frequencies with neutrals. The vertical 
perturbations along the background magnetic field have been investigated. We have included 
the effect of perturbation of coUisional frequencies due to density perturbations of species. 
We have shown that due to collisions of charged species with neutrals and neutrals with 
charged species the latter experience the back reaction on their perturbations. So far as 
the neutrals participate in the electromagnetic perturbations only because of collisions with 
the charged species, an adequate description of multicomponent plasmas including neutrals 
is to express the neutral dynamics through the dynamics of charged species and substitute 
induced velocities of neutrals into coUisional terms of the momentum equations for charged 
species. Then using Faraday's and Ampere's laws, we can derive the dispersion relation 
and/or investigate nonlinear structures. 
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The viscosity of magnetized species is important when k'^pj > 1, where pj — vrj/oJcj is 
the Larmor radius. The viscosity of neutrals is neghgible in comparison to the thermal pres- 
sure for the low frequency perturbations, uj <C Vnn- For the one- dimensional perturbations 
along the magnetic field the thermal pressure is present in the longitudinal perturbed veloc- 
ities of species. In the transverse velocity of neutrals, the viscosity of neutrals is important 
when Xnz ^ J^ni or /c^c^^ 3> I'ni^nn- Howcvcr, the neutrals are weakly coupled with ions in 
this case and are immobile in perturbations. 

The growth rates of perturbations in the wide region of the wave number have been 
found. We have derived that the long wavelength part of spectrum can be excited in medium, 
where CAi 3> Cgi 3> c^, and the short wavelength perturbations are excited at Cgi S> c^j, where 
Csi-, ca, and CAi are the ion thermal velocity, Alfven velocity including the mass density of 
neutrals, and ion Alfven velocity including the ion mass density, respectively. We have 
shown that the viscosity plays the role only for the short wavelength edge of spectrum when 
the neutrals have a weak collisional coui^ling with ions. In the cases of strong neutral- 
ion coupling, the viscosity is negligible. The resistivity is absent in the one-dimensional 
perturbations along the magnetic field. 

On the simple example, we have demonstrated that there is discrepancies between the 
standard MHD results and multicomponent approach at consideration of the resistivity. 
From the latter approach it is followed that the resistivity has an anisotropic nature. If the 
resistivity can play an important role in damping of perturbations, this effect can results in 
anisotropic turbulence. 

Electromagnetic streaming instabilities considered in the present paper can be a source 
of turbulence in thermal regions of astrophysical objects. 

The insightful and constructive comments and suggestions of the anonymous referee are 
gratefully acknowledged. 
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